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Abstract 

Naive dimensional analysis (NDA) is a widely used ansatz to estimate coupling con- 
stants among composite states emerging from dynamics of a strongly coupled gauge 
theory. However, the validity of NDA is still unclear because of the difficulty in calcu- 
lating these quantities in strongly coupled theories. We examine the NDA ansatz using 
gauge/string duality, by estimating glueball coupling constants from gravitational de- 
scription. The NDA scaling rule for coupling constants of some types of glueballs is 
verified and extended by both generic estimation and numerical calculations. The scal- 
ing rule verified in this article can be applied to some class of quiver gauge theories as 
well, not just to gauge theories with a single gauge group SU(N C ). 



1 Introduction 



It is beneficial to theoretically understand parameters in the low-energy effective theory of 
hadrons which emerges from a strong dynamics of a gauge theory. Models beyond the stan- 
dard model sometimes contain a strongly coupled sector (e.g. technicolor models, dynamical 
supersymmetry breaking models, etc.), and the parameters of the low-energy effective theory 
of hadrons can be observable or at least relevant to phenomenology. Even within the stan- 
dard model, the chiral Lagrangian is an effective theory of QCD. Although the parameters of 
the effective theories such as masses and coupling constants should be determined in terms 
of parameters in the theories at short distance, it is often difficult to calculate them due to 
strong coupling. 

However, there is an ansatz that is known to be reasonable to some extent; it is called 
naive dimensional analysis (NDA) [H 121 |3l H] , which roughly guesses magnitude of coupling 
constants among hadrons. In the NDA ansatz, the effective action of glueballs is given by 



where Anda is a parameter with mass dimension one and fields collectively denoted by <ft(x) 
represent glueballs. All terms in C((j)(x), <9 M , A ND a) are assumed to have dimensionless coeffi- 
cients of order unity (apart from an appropriate power of A ND a), and consequently, it follows 
that Anda is the mass scale of hadrons (except for Nambu-Goldstone bosons). The overall 
factor should be replaced with N c in the case of effective action of mesons (5j [6] , so that 
the N c scaling rule in a large N c gauge theory is satisfied. The essential point of the NDA 
ansatz is that the overall factor contains 4tt (~ 13), which is sizable compared with unity; 
this overall prefactor may not precisely be ^/(47r) 2 , but the NDA ansatz assumes that the 
deviation — represented by a yet undertermined factor 1 / /3 2 — is of order unity. After rescaling 
of 4>(x) so that the (f>{x) fields have canonically normalized kinetic terms, the coefficients of 
the interaction terms become of the forms, for example, 



for arbitrary integers I (> 1); any (/ + 2)-point coupling constants will scale in / as 
(47r/3/A c ) 7 x (9(1) in unit of Anda- There is an experimental support for the NDA ansatz (pQ 
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the coupling constant of three-pion-interaction (dir) 2 ir is f~ l in the chiral Lagrangian, and 
it is reasonably close to the prediction of the NDA ansatz, An 1 da ^-^72^, when experimental 
values, f n ~ 100 MeV , Anda ~ rn P ~ 1 GeV and N c = 3 are substituted. 

There are a couple of different ways to argue in favor of the NDA ansatz ([T]) also from 
theoretical perspectives. One of them is to require an ansatz of "loop saturation" in loop 
calculation in the effective theory. The loop saturation ansatz requires that each loop con- 
tribution to an amplitude is comparable to the tree level one, if the external momenta are 
taken to be of order Anda and the UV-cutoff of loop momenta is also set to be Anda- It 
follows from this ansatz that (/ + 2)-point coupling constants among arbitrary hadrons are 
0((47r) 7 ). However, the result seems to be incompatible with the N c scaling law, because 
the (/ + 2)-point coupling constants of glueballs should be proportional to N~ ! in a large N c 
gauge theory. If one multiplies (47r) 7 by adequate powers of N c in order to satisfy the correct 
N c scaling law, one obtains the NDA ansatz ([T]). 

Another argument for the ansatz (P, [2]) is based on the large N c expansion. Suppose that 
an operator O is a single trace operator such as Tr(FF) which can create glueball states \h n ), 
(vac|0(x)|/i n ,(p)) = f n e ipx . The (/ + 2)-point correlation function is given in planar limit by 

N 2 

(0( Pl )0(p 2 )...0(p I+2 )} = (2vr)V(5>)J /+2 (A;^)-^. (3) 

Z — / 1071" 

i 

If the 't Hooft coupling A = gy M N c is small, Ji + 2(\;Pi) can be calculated perturbatively in 
the gauge theory and it is given schematically by 

2 



Ji +2 (\; P i)~P d{I+2)+4 



0(l) + 0(-^-]+0 / 



167T 2 / \ 16n 2 



(4) 



where d is the dimension of the operator 0(p). The last factor Nl/(lQn) 2 comes from two 
color loops of gluons and a loop factor. For large A, perturbation breaks down and we do not 
know the form of Jj +2 (X]Pi). If we knew J/ +2 (A;pj) for large A, we could read mass spectra 
m n , decay constants /„, and coupling constants g ni ... ni+2 of glueballs from the behavior of 
Ji + 2(\]Pi) around mass poles p 2 — > —m^. of each external momenta; 

N 2 f 2 N 2 1+2 ( f \ 
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If we assume that the unknown functions Ji + 2(\',Pi) are of the form p d ( I + 2 )+ 4 x 0(1) as in 
the case of small 't Hooft coupling, we may obtain 



where m is the typical mass scale of glueballs. The second expression is in accordance with 
the NDA ansatz 

Both of the above two arguments, however, are far from being a clear justification for the 
NDA ansatz. The loop saturation in the first argument is an ansatz in itself. Moreover, the 
apparent contradiction between the loop saturation and the N c scaling may cast doubt on 
the reasoning of loop saturation. One may not be satisfied with the second argument either, 
because the estimation for the Green functions J/ + 2(A;p;) from perturbative calculation has 
no justification for a large 't Hooft coupling. 

In this article, we study the NDA ansatz (CQ) by means of gauge/string duality [7J El |9] for 
some types of scalar glueballs. Although the NDA has been checked only by experimental 
values of etc. in QCD so far, the gauge/string duality provides us with plenty of strongly 
coupled gauge theories which are calculable in the gravity duals. We leave a study of the 
NDA for other types of hadrons (such as glueballs with higher spin and flavor non-singlet 
mesons) for future work. The rest of the paper is organized as follows. In section 2, the 
glueball coupling constants are given in terms of the language of a string theory. They 
include overlap integrals which we cannot calculate analytically. We will verify the NDA by 
giving a generic estimation for the overlap integrals. In section 3, we confirm the estimation 
and the NDA scaling rule by numerical calculation in two specific gravity backgrounds. In 
section 4, we briefly summarize our study. 

Mass spectra and several three or four point couplings among composite states have 
already been studied by means of gauge/ string duality in [TTJ], [TTJ [12j [131 EJ EES HE] in several 
gravitational models. Our aim in this article is to find a rule governing the behavior of 
large numbers of coupling constants, not only three or four point couplings but generally 
(/ + 2)-point couplings. 

Note also that we will not only verify the NDA ansatz (pfl 12]), which is for SU(N C ) gauge 
theories. In this paper, we derive an extended version of NDA scaling rule, which cover more 
general gauge theories (such as a class of quiver gauge theories where the gauge group contains 
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p copies of SU(N C )). In the extended version of the NDA scaling rule, N c in the overall factor 
in (JTJ) will be replaced with y/a multiplied by an 0(1) factor, where a is something like a 
"degree of freedom" of a gauge theory. 

2 Glueball Coupling Constants in Gauge/String Dual- 
ity 

In the gauge/string duality [7J El [9], some gauge theories with a large number of color N c 
and a large 't Hooft coupling A = gY M N c are dual to string theories on warped spaces. A 
gravity dual of a four- dimensional conformal field theory is the type-IIB string theory on a 
spacetime AdS§ x W with a metric 

ds 2 = R 2 /z 2 (ri^ u dx^dx u + dz 2 ) + R 2 ds 2 w , < z < oo, (7) 

where W is a five- dimensional Einstein manifold. An infinite number of examples of W are 
known, with each manifold W corresponding to a conformal field theory. For example, the 
Type IIB string with W = S 5 is dual to M = 4 SU(N C ) super Yang-Mills theory [7], and 
a choice W = T 1,1 corresponds to an M = 1 SU(N C ) x SU(N C ) gauge theory with four 
chiral multiplets in the bifundamental representation [T7j. In references [HJ [19], one can 
find an infinite number of examples of gravity background geometries with W = Y p,q , which 
correspond to (SU(N C )) P = SU(N C ) x SU(N C ) x • • • x SU(N C ) (p factors of SU(N C )) quiver 
gauge theories with several chiral multiplets in the bifundamental representations between 
two of the p SU (N c ) groups. 

Gravity duals of confining gauge theories have also been constructed [201 [ZD [221 [231 121] • 
In most of them, the metric background in the string frame may be written without loss of 
generality as, 

ds 2 = {fiz^i^dx^dx" + dz 2 ) + R 2 ds 2 Wz , (8) 

with an appropriate definition of the ^-coordinate. The warp factor R 2 /z 2 in (JT]) for the case 
of conformal theories is replaced by a more general form (f(z))~ 2 . The internal manifold 
W z can also be dependent on z, and the other supergravity fields generally have nontrivial 
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background. Although details of the background are model dependent, most of the gravita- 
tional models which are dual to confining gauge theories have at least one thing in common; 
the radial coordinate terminates at finite valuqj z = z max - Normalizable wavef unctions of 
the fluctuations of supergravity fields around their backgrounds become discrete like cavity 
modes when the range of the radial coordinate (0 < z < z max ) is finite. Each composite state, 
i.e. glueball, in a gauge theory is described by each normalizable wavefunction in the dual 
gravity theory. The quantity z max corresponds to the mass scale of glueballs in the gauge 
theory. 

The four dimensional effective action of glueballs is obtained by dimensional reduction: 
expanding the fields in the ten dimensional type-IIB supergravity action in terms of the 
complete set of wavefunctions on the radial and compact directions and integrating over 
those directions. We can carry out this procedure by two steps of dimensional reduction; the 
first step is to construct a five dimensional supergravity action by integrating the complete 
set on the five dimensional compact manifold W z , and the second step is to obtain the 
four dimensional effective action of glueballs by integrating the complete set on the radial 
2-coordinate. 

2.1 An Intermediate Description in Five Dimensions 

In order to examine how glueball couplings are controlled in gauge/string duality, we start 
from dimensional reduction on five dimensional internal manifold in the case of the conformal 
background (JTj). It is useful to study the conformal case first, before going to the confining 
case (jHJ), because the conformal background ([7j) approximates the confining background (jSJ) 
to some extent, except for deep IR-region z ~ z m3iX . 

In the case of the conformal background AdS^ X W ([Tj), dilaton has constant background, 
= g s , and the self dual five form F 5 — (1 + *)JF§ has background T§ = ^Tpj^j ^ e vol(W), 
where vol(VT) is the volume form of W and Vol(VF) is its volume; Vol(W) = J w vo\(W). The 
backgrounds of the other fields, Co, B% and C2, are constant, which we assume to be zero 



lr The necessary and sufficient condition for confinement of heavy quarks are that (/(z)) -2 has a nonzero 
minimal value [25j . but this condition necessarily provides neither a finite z max nor discrete spectrum of 
glueballs. However, in most of gravitational models with confinement of heavy quarks, z max is also finite and 
(/(z)) -2 has its finite minimal value at z = z max . 
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for simplicity. After dimensional reduction on the five dimensional compact manifold W, the 
supergravity action using the Einstein frame metric becomes 



S = (fl5v °' ( r» x R3 [** r % ( + «w> - -yw + (ftcn ) + . . . 



^ldfis V Jo z 



(9) 



where 2k^ = (27r) 7 a /4 is the gravitational coupling constant of ten dimensional type IIB 
supergravity. Here we only keep track of 4>(x, z) and c(x, z) fields in five dimensions^, which 
correspond to fluctuations of $ (dilaton) and Cq (RR scalar) with a constant profile on W, 
respectively. The AdS radius R is given by 



>4 ..A.. a 1 ... .aVolfS 5 ) 



^ = te «« a Mf)= 4 *° Vo#f- < 10) 



Then the overall factor is rewritten as 



1 An 

x (i? 5 voi(H/)) x r 3 = -4, (n; 



with the definitions 



"SB- 

The supergravity action (jHj) becomes an expression that is interesting in the context of 
the NDA, by canonically rescaling the fields, 



S = oTx 



dz 



-Imr+id.m-liidcr+id.c') 2 ) 



(13) 



Here, <fr' = <p — (<&) and c' = g s c. The second line plays a role of interaction terms, which 
comes from the exponential dilaton factor e 2<? \ Besides combinatoric factors, the (J + 2)- 
point interaction terms of the five dimensional fields <p' and d have coefficients ((47r) / y/2a) T . 



2 The terms + • • • in Q represent terms including the fluctuations other than <f> and c. However, it 
includes no mixing in the bilinear terms between or c and the other fluctuations in the case of the AdS 
background considered here. On the other hand, the terms + • • • also includes interaction terms among <f> 
or c and the other fluctuations such as (d m (j))(d n (j))h mn , where h m „ is a fluctuation of the metric. We are 
keeping only the interaction terms among </> and c. 
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Although these coefficients are not glueballs coupling constant themselves, this expression 
f fT3|) is already quite indicative that the coupling constants of glueballs may show the behavior 
predicted by NDA (TJ. 

We saw that the parameter a in (I12p controls the interactions of the closed string fields. 
The parameter a corresponds to the central charge a of the conformal field theorjn which is 
dual to AdSs x W [27J 12H] • The central charge a represents something like degrees of freedom 
of the conformal field theoryfl For example, the property a oc N 2 can be understood because 
degrees of freedom of gluons and other fields in the adjoint representation of SU(N C ) gauge 
group are 0{N 2 ). Moreover, in the gravity dual of (SU(N C )) P quiver gauge theory, p' defined 
in (112j) is 0(p) [T8l [19], so we have a ~ pN 2 . This is also reasonable for the interpretation 
that a is roughly degrees of freedom of the conformal field theory because degrees of freedom 
of gauge theory with the gauge group SU (N) p should be about pN 2 . We also find that a has 
a lower bound 4a > N 2 [2H] because it is proven that p' has a lower bound p' > 1 in general 
Einstein manifold. Therefore, the parameter a makes the couplings in ffl3l) always smaller 
than (47r) 7 . 

2.2 Effective Action in Four Dimensions 

In order to obtain the effective action of glueballs on the four dimensional spacetime, we 
also have to integrate the fifth coordinate z. The supergravity action on AdS§ (0 [T3|) is 
no longer suitable for this purpose, because conformal theories do not give rise to a discrete 
hadron spectrum. We start to consider supergravity on confining geometries (JH1) instead, and 
describe coupling constants of glueballs. 

Let us first consider a theory which is asymptotically conformal in the UV region. In this 
case, only minor modification to the five- dimensional description (JH1 USD is necessary. In the 
gravity dual language, (f(z))~ 2 ~ (R/z) 2 in the z — > limit, and the central charge a in the 
UV limit is defined by ( IT2|) by using Vo\(W z ) at z = 0. On such a geometry, with the same 



3 For a review, see e.g. In conformal field theory, another central charge c is also defined. Although 

a and c are different in generic conformal field theories, one has the relation a = c when the theory has a 
weakly coupled gravity description as in the case we are studying in this paper. We just use the symbol a 
throughout this paper. 

4 For example, in a free massless field theory with n s scalars, rif Weyl fermions and n v vector bosons, one 
has a = (n s + ^-rif + 62n„) and c = —^j (n s + 3nf + 12n v ). 
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rescaling of the fields as in f )13p . the five dimensional supergravity action is written as 



lax d ~ 

S= I (fx I ^rY(z) 



z 3 



\m'f + {d-A'f)-\{{ddf + {d z c'f) 



2!/! 7 

where Y(z) is a dimensionless function defined as 



(14) 



y(z) - (/(^))" 3 (15 ) 

{ ) ~ {R/zf Vol(W^o)' (15) 
which is unity in the case of conformal geometries. In the z — > (UV) limit, Y(z) — >■ 1, and 
the integrand of ( |T4|) becomes identical to one of ( fi~3l) . 

In the rest of this article, we omit the terms in + ■ • • in (fl4|) . In general confining geometry, 
dilaton, RR scalar and the 3-form flux if 3 and F 3 would also have nontrivial background. In 
this case <f)'(x, z) and c'(x, z) would have mixings with other string fields such as the metric. 
We neglect this technical complexity in this article. We expect that such details of the IR 
background will affect coupling constants of glueballs at most by 0(1) factors, and hence 
they are unessential in trying to verify the NDA ansatz, whose predictions always come with 
uncertainty of order unity. 

We denote four dimensional glueball fields of the n-th excited modes by <f) n (x) and 
c n (x), which are created by operators dual to five-dimensional supergravity fields <f)'(x,z) 
and c'(x,z), respectively. We assign mass dimension one for glueball fields 4> n {x) and c n (x) 
just as usual for canonically normalized scalar fields in four dimensional field theory. The 
five dimensional fields (f)'(x,z) and c'(x,z) are decomposed into independent modes in four- 
dimensions, each one of which corresponds to a normalizable wavefunction ij) n (z)] 



(j)'(x, Z) = ^ i> m (z)<l) m (x), C'(X, Z) = ^ ^n{z)c n {x). (16) 
m=l n=l 

The normalizable modes ip n (z) are defined as the solutions of the eigen-equation given by 

z 3 Y~ 1 (z)d z (z- 3 Y(z)d z Mz)) = rr&Mz), (17) 
with the normalization condition 

lax 

—Y{z)ip n {z)il) m {z) = 5 nm , (18) 
z 6 
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which makes all the fields <p n (x) and c n (x) canonically normalized in the 4D effective action. 
The modes ip n (z) satisfies an appropriate IR-boundary condition which is imposed so that 
the field configuration in ten dimensional spacetime should be smooth at z = z max . The 
eigenvalue m n is the mass of glueballj^l 0n> c n . The effective action of glueballs is obtained 
by substituting (fT6D to 

The interaction part of the effective Lagrangian £ int includes (/ + 2)-point interaction 
terms which consist of two types of couplings; one is of the form 0^(9c n ) 2 , and the other 
r m Cn without derivatives: 



£* = -E E 

1=1 n\ti2m\. ..mi L 



A NDA 

a nW2,mi...mi 2\l\ ^ >mi ' ' ' < ^ mi 

(^C ni )(^ 



K 2 ' 1 ~ 

+b (I) NDA i 

1 "nin2,nii...mi 2\J\ 1 



(19) 



The coupling constants ai ltl2mi ... m , and % 1 n 2 ... mi ...m / are given by the following overlap integral 
of normalizable wavef unctions; 

ilm,...m, = ~7= X A NDA / ~T Y {z^m^n^mi ■ • ^m,, (20) 



nin 2 m 1 ...m I \ /q^ J N "NDA ^3 



/ 47T \ / (XX 

b { nL mi ... mi = [-7= x Kda / -^Y{z){d^ ni ){d z ^ n2 )^ mi . ..^ mi . (21) 



An \ r _ 2 /" 2max 

NDA 7o ^ 

These coupling constants have been made dimensionless by multiplying appropriate powers 
of a parameter Anda with mass dimension one. When we choose the mass scale Anda as the 
mass of the lightest glueball mass, these coupling constants are expected to be (D(A7rf3/N c ) 
by the NDA ansatz (El). 

We need to estimate the overlap integrals of fl20l I2T1) in order to examine whether there 
exists a rule governing the coupling constants of four dimensional effective theory like the 
NDA ansatz or not. The prefactor (An/y/^a) 1 is identical to the coefficient of (/ + 2)- 
point interaction terms of f JT3|) . and is just determined only from conformal region in UV, 
whereas the remaining overlap integrals are dependent on the detail of the IR geometry. 
The estimation of the overlap integrals unavoidably requires numerical calculations in each 



5 The normalizable wavefunctions ip n {z) and mass spectra m n are common in <fi' and c' in our study. This 
is because we are ignoring any kind of nontrivial vacuum expectation values and effects of mixing for the 
fields of supergravity, and therefore equations of motion for <p' and d become the same. 
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geometry and they will be the subject of section 3. Before numerical calculation, however, 
we will give a crude estimation independently of the detail of geometries. 

We may roughly estimate the overlap integrals for low excited modes and not too large I by 
approximating the integrand as a constant value in the IR region. On general UV-conformal 
geometry, normalizable wavefunctions behaves as z A (A is the conformal dimension and 
A = 4 for <j)(x,z) and c(x,z)) in the small z region, so the small z region have only small 
contribution to the overlap integrals. They oscillate with relatively large amplitudes in the 
IR region z ~ 0(z max ), say, z > z max /2. The normalizable wavefunctions of the first few 
excited modes have only a small number of nodes. Then it might be justified to approximate 
the integrand by a typical constant value in the IR region. Using the value of the integrand 
at around z ~ 7-2 ma x (7 ~ C(r)) as the typical value, the overlap integrals are estimated as 

I Q'm^mi ...rnj I ~ 

( "75= A NDA^max [ Y 'OjVV^n^m! ■ ■ ■ *Pmj] z ^ Zmax > ( 2 2) 



4 \ 1 

\b£ n2mi ...mr\ ~ ( "7= KdA^L [Y (z)(d^ ni )(d 2 ^ mi . . .^ mi ]_. • (23) 



W2mi...mi\ \ I In DA max L \~J \~z rni j \~z r rn« • • • rmi Jz~-yz max 

Applying the same approximation for bilinear terms, we also obtain 

ip n (z ~ 7z max ) ~ [z~l x Y(z ~ 72 max )] , d z ip n (z ~ 72 ma x) ~ mip n (z ~ 7^ max ). (24) 

We are assuming that Hi, rrij and I are not large, because the integrand would oscillate rapidly 
if they were large and the above estimation would break down. Here, 7z max is a typical value 
of z around which the integrand is peaked, and should be near the IR boundary, say, 7 ~ 1/2, 
and m is a typical mass of the low excited modes (say, m = mi). Substituting ff24|) to ff22l 
we obtain 



\^"n\nim\...mj I 



4 \ 1 1 

) A N DA^max (Y(z ~ 7^max)) _1/2 



2a J 



1 "nm2mi . 



.mj 1 



in 



A 



NDA 



2 



4"7T \ -1/2 
A N DA^max (Y (z ~ 7Z max )) 



-1 1 



2aJ 



(25) 
(26) 



The reason has the factor of (m/A^A) 2 is that f !23p contains two (d z ip n ys. 

With these crude approximations, one can find that there is a scaling rule for coupling 
constants like the NDA ansatz ([T|). When we choose Anda as the typical mass of glueballs, 

A N da = m, (27) 
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both of the (I + 2)-point coupling constants are estimated as 7-th powers of the same factor, 

1 



n\ri2m\ ...mj I 



n2mi ...mj I 



47T \ 



2a/ 

All \ 



~ T^max 

i 



)) 



-1/2 



2a J 



(mz n 



1 in the second line. The factor 



where we have also used an approximation Y(z ~ 7-z m ax) 
(mz max ) is expected to be 0(1) because l/z max corresponds to confinement scale of dual 
gauge theory, but this factor turns out to be slightly large in numerical calculations as we 
will see in the next section. Eq. ( 128]) with the choice of ( 127|) is just the same as what the 
NDA predicts fl]), with the identification of the NDA scaling factor 



Att/3 

~n7 



An 
\/2a 



(mz n 



(29) 



The NDA scaling rule shown in f[2"9"j) has been generalized from the original NDA ([T]) ; the 
factor N c / (3 in the original NDA rule ([1]) is generalized into \/2a/ ( m -2max)- The identification 
is natural in an SU(N C ) gauge theory because the factor ^/2a/(mz ma _ x ) is actually 0(N C ), 
which corresponds to the assumption of the original NDA that (3 is 0(1). However, Eq. (I29j) 
implies that (3 can take an arbitrarily small value, because Aa/N^ = p' can be arbitrarily 
large, for example, in a quiver gauge theory (SU(N C )) P with arbitrarily large p. On the other 
hand, we also find that f3 has an upper bound of 0(1) value because p' is bounded as p' > 1. 

So far, we have only focused on UV-conformal theories, but it is possible to extend the 
derivation above in order to cover theories with weakly running couplings even in the UV- 
limit. To do this, note that the ^-dependent function a(z) can be defined in gravity side even 
in non-conformal theorie EDI: 



a(z) 



(30) 



One can see that a(z) approaches the value a defined in (fT2l) when the geometry approaches 
AdSs x W. In a theory which has a weakly running coupling in UV, the function a(z) varies 



6 In a gauge theory which has both IR and UV fixed points, the central charge a (which is equal to c in 
supergravity approximation) has a relation auv > am (see |31j and references therein). The z-dependcnt 
function a(z) in (|30[) is defined so that it decreases monotonically as the holographic coordinate z is in- 
creased l30l. 
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only slowly in z, although it may vary rapidly in the deep IR region z ~ z max . Then the 
above NDA scaling rule (1291) is expected to be valid also in such theories by replacing a with 
(2, a typical value of a(z); 



Here, a = a(z ~ 7^ m ax) ; and z ~ 7-2 m ax is a typical value of z around which the overlap 
integral is dominated. The typical value a has an ambiguity which is at most 0(1) factor 
depending on the choice of 7 because a(z) is slowly varying, and such uncertainty is already 
expected in the estimation of this section. 

3 Numerical Check 

In the previous section, we have verified the NDA scaling rule by a crude estimation of the 
overlap integrals. However, it is difficult to estimate the error within that crude estimation. 
If the error in the estimation of (129]) were greater than the factor of order Air, that error would 
be too large to claim that the NDA really holds. In this section we want to confirm the NDA 
scaling rule ( 1291 l3Tj) by calculating coupling constants numerically, and checking that they 
are approximately given by the expected value (128]) within the range of 0(1) factors. We will 
perform numerical calculations in two specific gravity models. Our results are in considerably 
good agreement with (129 1 13TI) . 



First, we estimate the glueball couplings in "hard wall" model [32], which is regarded as the 

simplest toy model of IR-confining and UV-conformal gauge theories. The hard wall model 
simply introduces IR-boundary z = z max to the geometry ([7]) which is dual to a conformal 
field theory. Thus, the 5D action of this model is given by (fl4|) with Y(z) = 1. In this 
article, we choose Dirichlet boundary condition at the IR boundary on the wavefunctions, 
i.e., ^ n (z mffl ) = (just for simplicity). 

In the hard wall model, the mass spectra and normalizable wavefunctions are well-known. 




(31) 



3.1 Hard Wall Model 
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The wavefunction ip n {z) and mass eigenvalue m n of the n-th excited state are 

lf>n(z) = — '-— -j^tt—y , m n = fan/ZjpBi, (32) 
^max "2w2,nJ 

where J<i{%) is the Bessel function of the second order, and j2,n the n-th zero of ^(^O; 
j 2 ,i ~ 5.1, j2,2 — 8.4, j2,3 — 11- The dimensionless coefficients of the (/ + 2)-point couplings, 
an}n 2 m 1 ...mi and bn}n 2 m 1 ...m I , can be calculated by using these wavefunctions in (12"U | |2"T|) . 

We calculated these coefficients numerically, and obtained geometric means a[^ p and b[^ p 
of the ensemble of the (/ + 2)-point coupling constants involving lower excited statesj^l The 
result shows that the one of the typical values of the (/ + 2)-point coupling constants, a^, 
remains almost an /-independent constant of order unity (c ~ 0.3-0.4), after an appropriate 
scaling behavior is factored out: 

a ty P ~ CX ( 4=( A NDA-Zmax) X 0.8 ] . (33) 



2a 

See Table [TJ A similar result was also obtained for the other typical value b[y p of the (/ + 2)- 
point coupling constants when we take Anda = mi- 5.1/z max . 

This numerical result shows that the crude estimates in the previous section are quite 
accurate on average, and there is no question about the scaling behavior of the (/ + 2)-point 
coupling constants now. It is worthwhile to note, further, that the scaling factor in ( |29i) in 
the scaling rule we established is somewhat larger — by a factor of miz max ~ 5.1 or so in the 
hard wall model — than what was expected in the naive dimensional analysis!®! 

Although we have seen so far that the geometric means of (J + 2)-point coupling constants 
of hadrons, a[y P and b[y P , follow the scaling rule, such a rule will be of little value if the 
individual coupling constants a^ 7 1 ' ) „, 2mi ... m/ and bn} n2 m 1 ...m I take values wildly different from 
their averages. Our numerical study shows that the typical range of \an}n 2 m 1 ...m I \ is from 
a typ x [e - ' 7 ' 11 ' 1 ] to a[y p x [e +0 "i ( n"], with [e "'"' ] ~ 2. in the hard wall model. Similarly, it turns 
out that [e <J ini>] ~ 3. Thus, the coupling constants are typically within the factor of 2 ~ 3 



4yp, crjP Q , b[y' p and a[^ b , we used an ensemble of \an} n2mi ... mi \ and 
\bn}n 2 m l ...m I \ for a given value of /, with 1 < n\ < ri2 < mi < ••■ < mi < 3 and 1 < ni < ni < 
3, 1 < mi < • • • < mj < 3. 

8 The difference between (Air) and (4-7r) x 5.1 may yield different phenomenology, when the scaling rule is 
applied to physics beyond the Standard Model. 
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Table 1: Geometric means a\ yp and b[ yp and the standard deviations <r lna and <r ln6 calculated 
numerically in the hard wall model. The standard deviation of In |a^n 2mi ... mj |'s, that is, crj^, 
is presented in the 3rd column in the form of ln[ exp(<j^ a ) ], so that the range of typical 
values of \an}n2m 1 ...mi \ can be easily read out. 





"typ 


( 4vr \~ j 

, ANDA^max ' 0.8 

\^F2a J 


a {1) 

a lna 


b {1) 


( 4tt V' 

, A ND A^max ' 0.8 

KV2li J 




/ = 1 

1 = 2 
1 = 3 
7 = 4 


0.3 
0.3 
0.4 
0.4 


ln[2.] 
ln[2.] 
ln[2.] 
ln[2.] 




( mi > 
^A N da ) 

( mi > 
^A N da ) 

( mi > 
^A N da J 

( mi > 

^a ND a y 


x 0.4 

2 

x 0.3 
2 x 0.2 

2 

x 0.4 


ln[3.] 
ln[3.] 
ln[3.] 
ln[4.] 



from their geometric means. With the scaling factor (47r) x (ANDA^max) being much larger 
than the typical difference among the individual couplings [e CT i° ], we see that the scaling 
rule of the averaged value (133]) contains valuable information (albeit statistical) on individual 
(/ + 2)-point coupling constants. 

3.2 Klebanov-Strassler Metric 

One of the flaws of the hard wall model is that the IR-region of the geometry is very ad hoc; 
the IR-boundary z max is introduced by hand. Such a crude treatment is meant only to be a 
simplest toy model imaginable that implements the two essential ingredients of IR confining 
models: i) finite range of the holographic radius, 3 z mSuX > z > 0, and ii) existence of the 
minimal value of the warped factor f~ 2 (z). It is not meant at all to be a faithful (and hence 
stable) solution of the equation of motion of the Type IIB string theory. 

In a full solution of equations of motion of supergravity, however, the IR boundary z = 
2 max is not a singularity of the background geometry; the spacetime geometry is smooth in ten 
dimensions, and the internal geometry W z smoothly shrinks at z = 2 maX ; and we encounter 
a "boundary" of the geometry only after the description on 10 dimensions is reduced to that 
on the five dimensions. 

In the following, we construct a toy model using the Klebanov-Strassler metric so that the 
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model captures the above nature of faithful solutions to the equation of motions. With this 
toy model, where the IR geometry is treated in a more appropriate way than in the hard wall 
model, numerical calculation is carried out once again, in order to check the validity of the 
crude estimation of the overlap integration in section [21 and also to see how much individual 
(J + 2)-point coupling constants are different from their average. Because the Klebanov- 
Strassler metric does not asymptote to a pure AdS§ x W metric for some 5-dimensional 
manifold W in the UV region (z ~ 0), but maintains logarithmic running, our toy model 
using the Klebanov-Strassler metric also serves as an ideal test case of how to deal with the 
running a(z) function in the NDA scaling rule. 

The Klebanov-Strassler background is dual to an M = 1 SU(N) x SU(N+M) quiver gauge 
theory with chiral multiplets in the bifundamental representations [22]. This gauge theory 
experiences a cascade of Seiberg duality, SU(N) x SU(N + M) — »■ SU(N — M) x SU(N). 
In the IR, the duality cascade ends when the gauge group becomes SU (M) and confinement 
occurs. In the UV, the cascade of Seiberg duality continues unlimitedly, and it has weakly 
running RG flow in the UV; this theory is not UV conformal. 

Let us briefly review the Klebanov-Strassler metric, focusing only on the aspects that we 
need in the following. The background metric is given by [22] . 

ds 2 10 = h~ 1/2 (T)dx 2 + h 1/2 ( T )dsl (34) 
^^[{dr) 2 + (, 5 ) 2 ] + cosh 2 (I) [(, 3 ) 2 + (s 4 ) 2 ] 

+ sinh 2 (£)[(^ + (, 2 ) 2 ]], (35) 

where r (> 0) is the radial coordinate, g 1 , . . . ,g 5 are basis of 1-forms on the compact five 
dimensional space, and e 2 ^ 3 is a parameter with mass dimension —1. The two functions 
appearing in ( 134] [35]) are given by 



4/3 



, (sinh(2r) -2r) 1 /3 2 2 / 3 f°° , xcothx-1,. , .„ . „ 

K{t) = V ) 3 { / , h(r)=a— dx smh 2x - 2x ^ 36 

2 1 / d smh(r) 4 J T smh x 

where a = <i(g s Ma') 2 e~ 8 / 3 is a dimensionless constant determined by the dilaton vev, RR 

3-form flux M and the parameter e 2 / 3 . We can rewrite the metric 034|) into the form (jSJ), by 

defining the holographic coordinate z and the function f(z) in ([8]) as 

£ = -$n$P ^ = = A*) = Wr))^. (37) 
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In the UV limit (r — > oo, z — > 0), the asymptotic form of K{r) and /i(r) are given by 
K(t) -> 2 1 / 3 e" T / 3 and /i(r) ct2 5 / 3 3re 4r / 3 , respectively, and hence the new holographic 
coordinate z is related to the original one r by z — > a 1 / 2 e 2 / 3 2~ 5 / 3 3r 1 / 2 e~ r / 3 asymptotically. 
The asymptotic form of the metric [33] is given (in this coordinate z) byj^l 

ds 2 w -> ^(ln^A)) 172 ^^^ + ^ 2 ) + i^On^/z)) 1 / 2 ^,,, i? 2 = ^sMa 1 ^ (gg) 
£ 2 v 2 

where ds^i,! = |(g 5 ) 2 + hYli=i(9 1 ) 2 * s the metric of a five dimensional compact manifold, 
known as T 1,1 , which is topologically S 2 x S* 3 . Thus, the Klebanov-Strassler metric is ap- 
proximately that of AdSs x T 1,1 in the UV limit, except for the slowly varying logarithmic 
factors hi(z /z); z is a constant of order a 1//2 e 2//3 . 

In the IR limit (r — > 0), the functions in ( 1361) have finite limits, i^(0) = (2/3) ^ and 
h(0) ~ a x 0.285. The two directions g 1 and g 2 on the compact space shrink as oc \{dr) 2 + 
tK^ 1 ) 2 + (fi' 2 ) 2 ]- This form of the metric in the r ~ region is like a three dimensional flat 
metric written in polar coordinates, implying that the geometry ends smoothly at r = 0. The 
three other directions spanned by g%, and g§ form a three dimensional sphere at r = 0. 
Thus, the internal six- dimensional geometry is locally R 3 x S 3 , and r is a radial coordinate 
of R 3 [32]- In the z-coordinate f)37p . the smooth endpoint of the geometry r = corresponds 
to the maximal value of z, z max = z(r = 0) ~ (a 1 / 2 e 2 / 3 ) x 1.50. The warp factor (f(z))~ 2 has 
its nonzero minimal value xmxy[(f (z))~ 2 ] = (f(z ma _ x ))~ 2 — (ax 0.285) _1 / 2 at z = z max , which 
can also be written as min[/(z)]~ 2 ~ 2.65 x (-R/z max ) 2 - 

Let us now carry out dimensional reduction of Type IIB supergravity action on the 
Klebanov-Strassler metric (135]) . first to five dimensions. In non-conformal theories like this, 
the factor (R 5 Vol(W z )) x f~ 3 (z) most relevant to the dimensional reduction of dilaton is not 
simply proportional to 1/z 3 , and it does not even asymptotes 1/z 3 at small z in the UV 
non-conformal theories. In the Klebanov-Strassler metric, however, 



[(R 5 Vol(W z )) x f~ 3 (z 



/ox 5/2 

( I J ^(r)/i 5/4 (r) sinh 2 (r) x e 10/3 Vo\(T 1 ^ 



x Vol(T^) x ^ Z j z)) ' 

yd 




3 Wc follow the definition of R 2 in 
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in the UV region, and it is a reasonable approach to take only the i? 8 Vol(T 1,1 ) factor out 
of the integral, just like in (Q, while keeping (\n(zo/z)) 2 and 1/z 3 inside the integral. The 
action in five dimensions starts with the following terms: 

(40) 

4>'(x,z) and c'(x,z) are dilaton and RR 0-form fieldsj^j respectively, with a constant profile 
over the internal manifold W = T 1,1 . Here, 

f-\z) (^Vol(^)) 
1 ' (R/z) 3 R 5 Vo\(T^) ' 1 ' 

/»ra » T (3g s M 2 \ , Vo\(S 5 ) 27 



2tt y ' * Vol(Tbi) 16' 

One can see in a straightforward computation that the a(z) function defined in ( 130]) has 
an asymptotic form a(z) —> a c g (ln(z /^)) 2 in the UV (small z) region, precisely with the 
coefficient a c s defined above. 

In the rest of this article, we ignore all the terms denoted by + • • • , and take (1381) without 

the "H terms" as the starting point of a toy model. Non-trivial 3-form flux background H3 

and F 3 in the Klebanov-Strassler solution would potentially generate potential and mixing of 
4>' and d even at this level of description at five dimensions, and consequently would affect the 
mass spectra and coupling constants of hadrons in this theory. All these effects are ignored, 
however. Thus, numerical results from this toy model should not be taken literally as results 
of Klebanov-Strassler solution of string theory. We primarily use this toy modej^] for the 
purpose we already stated at the beginning of this section 13.21 

Now the effects of the confining geometry in the IR region (z ~ z max ) and the logarithmic 
violation of conformal symmetry even in UV region z — > are all encoded in the measure 
function Y(z) in this model. The profile of Y(z) in Figure [1] is neither constant (as in 
conformal theories) nor approaches a constant value in the UV (z — > 0) limit, reflecting the 



10 The field redefinition just below (|13[) does not involve any technical subtleties, because the dilaton vev 
is constant in the Klebanov-Strassler background. 

11 One might also expect that the numerically calculated coupling of hadrons of this toy model will provide 
a decent guess of those in the Klebanov-Strassler model within an error of a factor of order unity, which is 
good enough for the main subject of this article, but we prefer to make an error on a safe side, and would 
not push the argument that far in this article. 
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Figure 1: Y(z) in ( 140]) using the Klebanov-Strassler metric. 



nature of the Klebanov-Strassler metric fl35l) . Although the measure Y(z) of the overlap 
integration diverges as ()n(zo/z)) 2 in the UV region z — > 0, and vanishes as oc (z — z max ) 2 
at the IR boundary, its value remains quite stable and moderate, 4 ~ 1, for a middle range 
0.4 < z/z max < 0.8 in the holographic radius. This is where dominant contribution to the 
overlap integrals comes from, and hence there is no reason to suspect that the discussion in 
section [2] might go wrong in this case; in evaluating a(z) at z ~ 7z ma x to be used as a in the 
NDA scaling factor (I3ip . we only need to take 7 in this middle range, 0.4 ~ 0.8. 

Numerical calculation of hadron couplings in this model also shows that the (/ + 2)-point 
coupling constants follow the scaling rule, when the scaling factor is chosen to be 



See Table El Because a evaluated somewhere in the middle range 7 G [0.4, 0.8] of the holo- 
graphic radius is much the same as a c g numerically, the crude (and model independent) 
estimation of the overlap integral in section [2] holds true quite accurately for this model. 

It is also worthwhile to note, just like in section 13.14 that the factor (ANDA^max) becomes 
6.2, if the lightest hadron mass eigenvalue mi is used for Anda in this toy model. The NDA 
scaling rule holds, but the scaling factor, (An) x 9 apart from the degree of freedom factor 




(43) 




, is larger than expected in the NDA ansatz. 
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Table 2: Geometric means a[y P and 6^ and the standard deviations of ln|ai^ 
di) i Ji) n „ A JX) ' 



and 



\n2ra\...mj \ 

In b 

I' a ina and °"in6: respectively, calculated numerically in a model (J3DJ) using the 



Klebanov-Strassler metric. Just like in Tabled! in the 3rd and 5th columns are presented 
in the form of ln[ exp(a^) ]. 



l typ 



47T 



ANDA^max ' 1-5 



(T 



(0 

In a 



fc (0 f _J^_ 
typ Vv^ 



max x • <-> 



a 



(0 
In 6 



/ = 1 

7 = 2 

1 = 3 

7 = 4 



0.3 
0.3 
0.4 
0.5 



ln[2.] 
ln[l.] 
ln[2.] 
ln[2.] 



mi 



mi 



mi 



mi 



x 0.3 



x 0.2 



x 0.2 



x 0.2 



ln[3.] 
ln[4.] 
ln[4.] 
ln[3.] 



4 Summary 

Our study based on gauge/string duality is summarized as follows: the NDA scaling rule 
does exist for the scalar glueballs dual to <fi and c fields, and the scaling factor is given by 
( 13~TT) . The error of the scaling factor is within a factor two or so (0.8 for hard wall model and 
1.5 for the model using the Klebanov-Strassler metric.) The uncertainty in the choice of a e s 
is also of the same order. The N c dependence of hadron couplings following from the large 
N c argument (and from string theory) is now generalized in the language of a(z) function 
defined in holographic models, which roughly characterizes the degree of freedom in the dual 
theories. With this generalization, the NDA scaling rule can be applied also to some class of 
quiver gauge theories. 

If we are to use the mass of the lightest (non-Nambu-Goldstone boson) hadron as the 
Anda mass scale in (JT]), the scaling factor comes out to be larger than the conventional one, 
by (mz max ), which is about 5 ~ 6 in the two models we studied numerically. We can say 
that the (dimensionless) coefficients of hadron couplings are systematically larger than the 
conventional NDA ansatz, but alternatively, we can also say, by taking Anda — l/^max; that 
the coefficients of hadron couplings are just as expected in the conventional NDA ansatz, and 
mass parameters are somewhat larger than naively expected from Anda- We have nothing 
to say, however, about whether the large value of mz mas is a generic feature of hadrons from 
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strongly coupled theories, or just happens to be a specific feature of the two models we 
studied. 

So far, we are not sure whether coupling constants of any other hadrons have the NDA 
scaling rule, and if so, what are the NDA scaling factors. These are subjects in our future 
work. 
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